Multiple dynamic transitions in nonequilibrium work fluctuations 
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The time-dependent work probability distribution function P{W) is investigated analytically for 
a diffusing particle trapped by an anisotropic harmonic potential and driven by a nonconservative 
drift force in two dimensions. We find that the exponential tail shape of P{W) characterizing rare- 
event probabilities undergoes a sequence of dynamic transitions in time. These remarkable locking- 
unlocking type transitions result from an intricate interplay between a rotational mode induced 
by the nonconservative force and an anisotropic decaying mode due to the conservative attractive 
force. We expect that most of high-dimensional dynamical systems should exhibit similar multiple 
dynamic transitions. 
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Systems in thermal equilibrium (EQ) are governed by 
the principle of statistical mechanics. It serves as the 
unified framework for the study of thermodynamic prop- 
erties of EQ systems, and has been successful since it was 
established centuries ago. On the contrary, such a prin- 
ciple, except for the second law of thermodynamics, was 
absent for nonequilibrium (NEQ) systems, which made 
it difficult to understand NEQ phenomena. Recently, 
discovery of the fluctuation theorem (FT) opened a new 
perspective on NEQ processes and has attracted a lot of 
interests from researchers. The FT refers to identity rela- 
tions for a thermodynamic quantity, such as work, heat, 
or entropy production, that are derived theoretically for 
a wide class of NEQ process [H-Q , and are also confirmed 
experimentally |10l4l2| . It not only serves as a criterion 
diagnostic to NEQ, but also sheds light on quantitative 
understanding of NEQ fluctuations |13l4l7 1 . 

Consider a system which is in thermal EQ with a heat 
reservoir at the inverse temperature (3. It is driven into 
a NEQ state by applying a nonconservative force [l^ 
or making contact to other reservoirs at different tem- 
peratures [H-dll. The NEQ work W done on the 
system over a finite time interval is known to satisfy 
the integral FT, known as the Jarzynski equality (JE) 
(e-P^) = 1 0,11, [ill . The JE, being combined with the 
Jensen's inequality, yields (W) > 0, which corresponds to 
the Kelvin statement of the thermodynamic second law. 
Furthermore, Crooks found that its probability distribu- 
tion function (PDF) P(M^)satisfies the detailed FT such 
as P{W)/P{—W) = [4|. It may deliver a surprising 
message that negative-VF events are possible, however 
they are exponentially rare and thus one still finds the 
average work (W) > 0. 

The FT evokes the importance of studying NEQ fluc- 
tuations, especially in the rare-event region (e.g. negative 
W) to measure the violation of the thermodynamic sec- 
ond law in a small system. Many studies have been done 
for the PDF's of the work and heat associated with NEQ 



jrocesses, theoretically and experimentally [l-12, 18- 
23j . These are to confirm the FT in the first 
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place, and then to investigate the nature of NEQ fluctu- 
ations. Normally it is a formidable task to find the PDF 
analytically for a specific NEQ process, so most studies 
are limited to special cases such as the large deviation 
study in the infinite-time limit (23 - [28} . 

In this Letter, we investigate the PDF of a NEQ work 
W over a finite time interval t in a two-dimensional lin- 
ear diffusion system (LDS), driven by a nonconservative 
linear drift force. The LDS is often referred to as a mul- 
tivariate Ornstein-Uhlenbeck process describing the mo- 
tion of a Brownian particle or a particle trapped by a 
linear force in the overdamped limit (iol . [sof . This sys- 
tem is simple enough that the PDF PjW) is analyti- 
cally tractable for finite time interval [ll,[3^ and still has 
many interesting experimental applications, e.g. a nano 
heat engine in contact with multiple heat reservoirs [l3| 
and a colloidal particle driven along periodic potential 
imposed by laser traps [isj . Our study reveals that even 
such a simple system displays a surprisingly rich dynamic 
behavior with a sequence of dynamic transitions in time. 

We briefly summarize the results with dimensionless 
W in units of 1//3: (i) The PDF has exponential tails 
with power-law prefactors as P{W) ~ g-w/w± 
in the W — > ±oo limit. The power-law exponents are the 
same in both sides (a-^. = a_ = a), and the characteristic 
works W+>0 and < satisfy 1/W+ + 1/W- = -1, 
which are required for the detailed FT. (ii) The power- 
law exponent a can take three different values of 0, 1/2, 
and 2. Accordingly, the PDF is categorized into type 
with a = 0, type I with a = 1/2, and type II with a = 2. 
Interestingly, W± continuously varies with time t for type 
and type I, while they are constants of time for type 
II. (iii) Typically, the system undergoes a dynamic tran- 
sition from type I to type II as t increases. The charac- 
teristic work W+ increases smoothly with t and suddenly 
becomes frozen at the transition time and afterwards. 
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This is a kind of a locking transition. More remarkably, in 
some parameter space, the PDF alternates between type 
I and type II indefinitely, i.e. infinite number of locking- 
unlocking type transitions. In general, a finite sequence 
of dynamic transitions is also possible as well as no tran- 
sition with type I in all time. Type is found without 
any dynamic transition, only in a special case. 
We consider a LDS with the equations of motion 



^9 c , t 



(1) 



where q = {qi, - ■ ■ ,qd)'^ is a d-dimensional vector, F = 
(Fij) is a constant positive-definite {d x d) force matrix, 
and ^{t) = (^i(t), • • • ,£,d{t))'^ is the Langevin noise sat- 
isfying 



(e.w> = o, (6(t)o(t)> = 2A,<5(t-t') 



(2) 



with a noise correlation matrix D = {Dij) which is sym- 
metric and positive-definite. After a similarity and a 
scale transformation, one can take the noise matrix as 
the identity matrix I (D = I) without loss of generality. 

One can decompose the force matrix into the conser- 
vative and nonconservative parts as F = Fc -|- fnc with 
Fc = F^. Fnc is not symmetric to generate a nonconser- 
vative force. When F„c = 0, the total force / = — F • q is 
conservative and can be represented as a minus gradient 
of a potential function V{q) ~ ^q^ ■ F ■ q. The steady- 
state PDF is given by the EQ Boltzmann distribution 
Peqig) (X e~^('^' satisfying the detailed balance condi- 
tion. A nonsymmetric force matrix (F„c 7^ 0) indicates 
the presence of a nonconservative force fnc = —^nc • Q 
which cannot be written as a gradient function. It drives 
the system out of EQ. Here, we only consider the anti- 
symmetric Fnc (= —^nc) foi' simplicity [3l[. 

Suppose that the system is prepared in the thermal 
EQ with Peq{q) oc e~(^/^)'i ''""''^ with the conservative 
force fc = — Fc • q only, then turn on the nonconservative 
force fnc at i = 0. The NEQ work (done by fnc) on the 
particle following a path ^(t) for < r < i is given by 



W[q(T)] = - f dT 
Jo 



dq{rf 

dT 



■Fr. 



(3) 



Defining the PDF P{W) = {5{W - H'[q(r)])) and its 
characteristic function g{t;X) = f dWe'^"^ P{W), we 
obtain ^ 

InGit; \) = - f dT Tr(A(r; A) - F) - ^ In 







det Fr 



(4) 

Here, the symmetric {d x d) matrix A is the solution of 

-^A(t; A) = -2A2 + AF + F'^A + A (5) 
dT 

with the initial condition A(0) ~ Fc, and the auxiliary 
matrices are defined as F(A) = F — 2AF„c and A(A) = 
(FJ^F-F^F)/2 0. 



We apply the formalism to a two-dimensional system 
with the force matrices, parameterized as 



l + u 

V 



1-u 
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with li > and £ > 0. One can set the trace of F to 
any positive number by the global rescaling of q and t. 
Here, it is set to be 2. Positive-definitencss of F and Fc 



requires that cr^ 



■v^ < 1. The equations of motion 



are covariant under rotation, and e and cr are invariant. 
So we may take v = without loss of generality. Note 
that the parameter e corresponds to the strength of NEQ 
driving (torque), while the parameter a represents the 
anisotropy of the harmonic potential. 

Before considering the general case, we present the 
solution in the special isotropic case with u = v = 0. 
In this case, the matrix A is proportional to the iden- 
tity matrix I as A = z{t;X) I. Then, Eq. ([S]) becomes 
dz/dT = -2z2 + 2z + 2e'^X{l - A) with z(0) = I and 
Eq. (|4]) is written as 



ln^;(t;A) = -2 / {z{t; X) - l)dT - In z{t; X) . (7) 
Jo 

The solution of the differential equation is given by 



A tanh( V At) 



anh(\/Ar) 



where 



1-Ae'^XiX-l) >0 



(8) 



(9) 



At A = 0, we have z{t;X) = ^ + ■2{t+i} • negative 
A, the solution in Eq. (|8]) should be continued analyt- 
ically using tanh(ie) = itan(e). Examining the above 
equations, we draw the following conclusions: (i) Since 
A(A) = A(l - A), one finds g(t] A) = g(t; 1 - A) at all 
t. This verifies the Crooks detailed FT. (ii) Given t, 
there exists Xo{t) > 1 such that z{t; X — Aq) ~ z{t; X ~ 
1 — Aq) = 0. The logarithmic term in Eq. ([7]) tells 
us that the characteristic function has a divergence as 
Q{t; A) ~ (Aq — A)^^ as a Aq from below. There is the 
same divergence at A = 1 — Aq. (iii) The simple poles 
manifest the exponential tails of P{W): 



P{W) 



e-(Mt)-i)w ^ ^ 00 
f,-Xo(t)\W\ ^ W -00 



(10) 



where Xo{t) monotonically decreases in t and asymptot- 
ically approaches to Ac = 5 + vT+e^/(2e) > 1. This 
PDF belongs to type 0. (iv) The large deviation func- 
tion e(A) = — limt^-cx) jlnC/(<;A) exists and is given by 
e(A) = ^A(A) — 1 for A(A) > and —00, otherwise. 

The analysis above gives us a lesson that the singu- 
larity of g, hence the tail behavior of P{W), is deter- 
mined from the root of det A{t; A) — 0. In the special 
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case {u = V = 0), two eigenvalues of A are degenerate. 
So, the root contributes to a simple pole of G and the 
pure exponential tail of P{W). When the degeneracy is 
broken, which is the case for any nonzero value of m or w, 
G has a square-root singularity at the root and one may 
expect the PDF of type I. However, the actual behav- 
ior is much richer and can be understood from the full 
solution of Eq. 

Here, we sketch briefly the way to find the exact solu- 
tion of Eq. ([5]). Details will be published elsewhere js^ . 
First, one needs to find a fixed-point solution As(A) such 



that 



rfA 
(It 



A=A^ 



0. It can be expressed as 



VF), 



(11) 



where an orthogonal matrix V (VV"^ = I) is to satisfy 
F^V^ - VF = F - F"^ to ensure = A^. Then, the 
matrix R(t; A) = (A — Ag)^^ satisfies a linear differential 
equation 



— = 21 - RFT - FR , 
dr 



(12) 



where F(A) = F — 2As and the initial condition is given 
by R(0) = (A(0) - As)"^ Its solution is 

R(r; A) = e-^^ R(0; A) e~^^' + 2 / dT'e'^'^er^'^^ 



= Roc +e-"^ (R(0;A)~Roo) e-"^^, (13) 

where R^o = (I + Q) with the antisymmetric matrix 
Q satisfying FQ + QF^ = F- F'^ [l^. Finally, one obtains 
the full time-dependent solution 



A(t;A) = A, + R(t;A)- 



(14) 



of Eq. ([5]) in any d dimension. Associated algebra is 
rather simple in two dimensions. We find Q = (F — 
F-^)/TrF and V is also easily determined. The explicit 
expression for Eq. is still rather complex, which will 
appear in 
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The solution is just the starting point for further anal- 
ysis of the PDF. We found that det A{t; A) exhibits a 
complex behavior depending on the parameter values of 
e, u, and v. There are four distinct cases, which are 
shown in Fig. [TJ The plots are obtained for a few values 
of u to a given e = 3/10. We took v = since only 
a ~ y/u^ + v^ is the relevant parameter. 

(i) When u is small enough (see Fig. [I{a)), the curve 
is not tangential to the r axis at any value of A. For 
any given t = t, one can find Aq such that det A(t, Ao) = 
(non-degenerate root). Then, detA(<, A) ~ (Ao — A) 
near A ~ Aq and, from Eq. (g]), G{t,\) - (Aq - A)-i/2. 
Hence, the PDF has a tail P{W) - \W\~^/'^e^<'^ in the 
W — >■ — 03 limit (type I), with Aq ~ Ao(t) monotonically 
decreasing with t to an asymptotic value Ac(u,e) > 1. 
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FIG. 1. (Color onlme) Time evolution of detA(r;A) when 
£ = 3/10 and (a) u = 0.65, (b) 0.7, (c) u* = ^109/200, and 
(d) 0.8 with several values of A > 1 (lower curves represent 
larger values of A). The curves tangential to the r axis (dot- 
ted horizontal line) are drawn with thick lines (red). The 
characteristic function Q is well defined only when det A > 0. 
Curves in the unphysical regions with det A < are plotted 
with dashed lines. 



(u) In the intermediate values of u (see Fig. [1Kb)), the 
curve is tangential to the t axis at multiple values of A. 
We will denote the time for the n-th tangential point as 
(marked with vertical solid lines) and the correspond- 
ing A value as A„. The curve, det A(t, A„), that is tangen- 
tial at T = i+ may cross the r axis at a later time denoted 
as r = t~ (marked with vertical dashed lines). This 
crossing is linear (non-degenerate root) and never hap- 
pens again later. Within the time interval t+ < t < t~, 
the characteristic function G{t;X) is finite for A < A„, 
and then it diverges discontinuously at A = A„. Hence 



the PDF has a tail P{W) ~ \W\-^e^"^ in the W ^ -oo 



limit (type II) [33|. Note that A„ is a constant of time 
within the finite time interval. Outside the interval, the 
PDF belongs to type I. Hence, the PDF alternates be- 
tween type I and type II many times as t increases. 

(iii) At the special value of m = m* = ^/{T+e^)/2, the 
curve, detA(T, A*), with A* = (1 + e)/{2e) is tangential 
to the r axis infinitely many times as 



detA(T;A*) 



e(l - £)(1 + e){l + coswt) 



4 e2^^ 



1+e 



COS UJT 



(15) 



with uj = \/2{l — e'^) and the tangential points at i+ ~ 
{2n - l)n/uj (?i = 1,2,.. .). Note that <,7 = t+^^. Hence, 
the PDF belongs to type I when t < (marked with 
the vertical solid line) and changes to type II afterwards 
except periodic instantaneous moments at t = {n > 
2) (marked with vertical dot-dashed lines). 

(iv) When u > u*, the curve is tangent to the t axis 
at a single value of A at r = ij*" (marked with the vertical 
solid line) without crossing the r axis later. Hence, the 
PDF belongs to type I for r < tf and type II afterwards 
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FIG. 2. (Color online) Constant-a; curves in the \o-t plane at 
e — 3/10. The value of u ranges from 0.55 (top) to 0.76 (bot- 
tom). The solid curves correspond to the PDF of type I, while 
the dashed line corresponds to the PDF of type II with the 
boundary drawn with the thick zigzag (blue) curve. The thick 
(red) curve corresponds to the case with u — u* . 
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FIG. 3. Phase diagram in the u-t plane for several values of 
e = 0.1 (bottom), 0.3, 0.5, 0.7, and 0.9 (top). The PDF be- 
longs to type II in the region surrounded by the phase bound- 
ary curve, and type I elsewhere except the line u = where 
the type-0 PDF is found. For £ > 1, the phase boundary 
disappears completely and no type II exists. 



forever. 

Our analysis reveals that the system undergoes a dy- 
namic transition in the tail shape of the PDF between 
type I characterized by P{W) ^ \W\-^^^e^<>'^ and the 
type II characterized by P{W) ^ \W\~^e^"^ for large 
negative W. The same applies for large positive W due 
to the symmetry of the detailed FT, see Eq. ((TO)) . The pa- 
rameter Ao , which corresponds to the inverse of the char- 
acteristic fluctuation size \ W- \ for the negative PDF tail, 
decreases continuously in time for type I, while it is locked 
for type II. For example, in case of (ii), Aq decreases up 
to t = , and is locked for a while till t = , then is 
unlocked a,i t = and decreases again till t ^ t2 , and 
so on. These locking- unlocking transitions occur many 
times as t goes by. 

The transition moments can be found numerically ex- 
actly from the roots of det A(t; A) = 0, using the solution 
given in Eq. (jl4p . We present the constant-u curves in 
the Xo-t plane at e = 3/10 in Fig. H The PDF un- 
dergoes a single dynamic transition for large values of 
M, infinitely many transitions for intermediate values of 
u, and no transition for small u. From the constant-w 
curves, one can construct the phase diagram in the u-t 
plane. The phase diagram is drawn in Fig. [31 

In summary, our analytic result shows that such a sim- 
ple linear diffusion system exhibits surprisingly complex 
NEQ fluctuations. It raises interesting questions for the 
mechanism of the dynamic transition. The conservative 
part of the drift force generates an anisotropic harmonic 
potential, which attracts the particle toward the origin. 
The nonconservative part acts like a torque, which drives 
the particle into a rotational motion. It is useful to intro- 
duce the polar coordinate to focus on the rotational dy- 
namics separately. Then, the dynamics of the polar angle 



may be written down effectively as d(j)/dt = a sin </)+fe+^, 
where the potential amplitude a should be proportional 
to the anisotropy cr and the constant driving force b 
should be proportional to the strength of the nonconser- 
vative force e. This is the equation of motion in a tilted 
periodic potential 15|. For a <^h (small cr), the particle 
has no time to relax in the potential well and drifts down 
into the steady state with a constant velocity. So there 
is no extra time scale except for the relaxation into the 
steady state. Thus, we expect no dynamic transitions 
and typc-I PDF forever. For a ^ b (large cr), the parti- 
cle sits at the potential well long enough and fully relaxes 
inside the well. Then, it can hop to the neighboring well 
due to the noise after a finite time, which can be deter- 
mined by the noise strength. This additional time scale 
exists and may set the transition time . As the parti- 
cle relaxation in the first well is fully developed already, 
there will be no additional time scale needed for succes- 
sive hoppings. Thus, we expect one dynamic transition 
from type-I to type-II PDF. When a w 6 (intermediate 
cr), it needs additional time scales for successive hoppings 
with incomplete relaxations within each potential well, 
which leads to multiple dynamic transitions. The above 
argument provides a plausible understanding of existence 
of multiple time scales, but does not fully capture under- 
lying mechanisms of dynamic locking-unlocking transi- 
tions. Since the radial component also fluctuates in our 
model, there are more possible complex routes to relax 
into the steady state. It is remarkable to find no smear- 
ing out of sharp dynamic transitions with fully locking 
states. We leave full intuitive physical understanding of 
these dynamic transitions for future works. 

As can be seen in our analysis, the remarkable char- 
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acteristics of multiple dynamic transitions should not be 
a pathological property of some special systems. Hence, 
we expect that any high-dimensional dynamical system 
driven by a nonconservative force should exhibit similar 
or more complex multiple dynamic transitions in a rea- 
sonably large parameter space. It would be interesting to 
observe these locking- unlocking features in experimental 
setups such as in 



19j by measurements in transient 
regime. 
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